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ABSTRACT 

A formulation to include the effects of wall-proximity in a 
second moment closure model is presented that utilizes a tensor 
representation for the redistribution term in the Reynolds stress 
equations. The wall-proximity effects are modeled through an el- 
liptic relaxation process of the tensor expansion coefficients that 
properly accounts for both correlation length and time scales as 
the wall is approached. DNS data and Reynolds stress solutions 
using a full differential approach at channel Reynolds number of 
590 are compared to the new model. 

INTRODUCTION 

The theoretical development, of higher-order closure 
models, such as Reynolds stress models, have primarily 
been formulated based on high Reynolds number assump- 
tions. Accounting for the influence of solid boundaries on 
these closure models has usually evolved from either a wall 
function approach or a modification to the high Reynolds 
number form of the pressure-related correlations and tensor 
dissipation rate, and predicated on the near-wall asymp- 
totic behavior of the various velocity second moments (e.g. 
Hanjalic 1994). 

A more broad based attempt to account for the prox- 
imity of a solid boundary is the elliptic relaxation approach 
introduced about a decade ago (Durbin, 1991), and fur- 
ther developed for second-moment closures (Durbin, 1993; 
Wizrnan, Laurence, Kanniche, Durbin & Demuren, 1996; 
Manceau and Hanjalic, 2000, Manceau, Carlson & Gatski, 
submitted for publication). In its two-equation form, the 
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v 2 — f model has been applied to a variety of flows (e.g. 
Durbin, 1995; Pettersson Reif, Durbin & Ooi, 1999). The 
new approach outlined here introduces a tensor function 
representing the combined effects of a near-wall velocity- 
pressure gradient correlation and anisotropic dissipation 
rate that asymptotes to a high Reynolds number form away 
from solid boundaries through an elliptic relaxation equa- 
tion for polynomial expansion coefficients. 


With the development of a generalized methodology 
for constructing tensor representations for the implicit al- 
gebraic turbulent anisotropy equations (Gatski fc Jongen, 
2000), it is possible to similarly construct tensor repre- 
sentations of only the redistributive part of the Reynolds 
stress equations (i.e. pressure-strain rate correlation and 
anisotropic dissipation rate equations). The elliptic relax- 
ation formulation outlined here can be viewed as an inter- 
mediate step between a fully explicit algebraic Reynolds 
stress model and a full differential Reynolds stress model. 
The intermediate step being the relaxation of the expansion 
coefficients of the tensor representation rather than that of 
the redistribution function itself. 


Comparisons are made with direct numerical simula- 
tion channel flow at a Reynolds number based on friction 
velocity of 590. The comparisons include the mean flow, 
turbulent kinetic energy and the Reynolds stresses for solu- 
tions from both the relaxation of the redistribution terms 
and relaxation of the corresponding expansion coefficients. 


Address all correspondence to this author. 
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THEORETICAL BACKGROUND AND DEVELOPMENT 

In this section, a mathematical framework is outlined 
that allows for the construction of a set of elliptic relax- 
ation equations for the polynomial expansion coefficients in 
a tensor representation of the redistribution terms in the 
Reynolds stress transport equation. 

The starting point in developing a elliptic relaxation 
formulation is the transport equation for the Reynolds 
stresses which can be written as 
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where the turbulent transport term, which included expres- 
sions for the triple-moment and the pressure transport, is 
modeled using a gradient-diffusion model ( C ^ = 0.15, and r c 
is a composite time scale defined in Eq. (9)). An anisotropic 
dissipation rate model has been introduced (e.g. Manceau, 
2000) into the redistributive term, so that sK fij is 
given by 


where a\ = 0.09. Note that this form of the dissipation 
rate equation (Durbin 1993, Manceau & Hanjalic, 2000) has 
introduced a composite time scale, r c , that replaces the ratio 
Kje in both the production and destruction of dissipation 
terms. 


Elliptic Relaxation Methodology: r — f Model 

The rescaled elliptic relaxation equation is driven by 
the high-Reynolds number form of the pressure-strain cor- 
relation, II, and a contribution from the Reynolds stress 
anisotropy, 2 £ c bij (away from the wall the dissipation rate 
is assumed to be isotropic dij = 0) 

(1 - L 2 V 2 ) fij = T (IT*. + 2e c b tj ) (6) 

where 



with the relaxation scales defined as 


eKfij = <j>ij - 2<r (dij - bij) (2) 


where the Reynolds stress anisotropy bij and the dissipation 
rate anisotropy dij are 
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The wall proximity effect is introduced through an elliptic 
relaxation of the redistributive model fij . 

Solutions are obtained for the Tn and T 22 normal 
Reynolds stresses, and the r \2 shear stress. The transport 
equation for the turbulent kinetic energy is used in lieu of 
the third normal stress 733 
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and the associated modeled transport equation for the tur- 
bulent dissipation rate, £, is 
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with <tk — 1-0, (T e — 1.3, C e 2 = 1.83, = 0.15 and 
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the present study and were adopted from previous elliptic 
relaxation formulations. 


In the current formulation, the goal is to implement a 
tensor representation for the relaxed redistribution function 
fij. As will be further discussed in the next section, such a 
representation would be consistent with a linear pressure- 
strain rate model. Previous implementations of the elliptic 
relaxation procedure (Manceau & Hanjalic, 2000) using the 
SSG pressure strain rate model (Speziale et al., 1991) uti- 
lized the full nonlinear form. The linear form of the SSG 
model implemented here is given by 


n h = -(c?£ c + C\v) b + KCnS 
+/\C 3 (bS + Sb - j{bS}l) - KC4 (bW - Wb) 

( 10 ) 


c;, = 1.44 ( 1 + «i T 


with Co 1 = 1.8, Cl = 3.4, Cl = 0.37, C 3 = 1.25, C A = 0.4. 
Since this linear form is used, the length scale constant Ci 
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needed to be modified from the previously used value. 

In the vicinity of the wall, the redistributive term is 
balanced by the viscous diffusion of the Reynolds stresses 
resulting in the boundary conditions for the fij (see Ta- 
ble 1). Only the 22- and 12-components of f have determi- 
nate solutions to the near-wall balance of the stress trans- 
port equations. Manceau et al. (submitted for publication) 
have used /n = fa = —fa/ 2 as boundary conditions on 
the remaining components to ensure that f is traceless. 

Table 1. RELAXATION TENSOR BOUNDARY CONDITIONS. 


Component 

Wall 
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“ 20,12 1(1) 

ely 4 {i) 

~ -zf22,w 
-20 ut?| (1) 

ely 4 (1) 

symmetry 

fa 

symmetry 

fa 

symmetry 

/l2 

0 


Representations and Elliptic Relaxation: t — j3 n Model 

While the elliptic relaxation formulation has already 
been applied within a full differential Reynolds stress model, 
a question arises about the role tensor representations can 
play within the framework of the elliptic relaxation proce- 
dure. 

As a first step, consider only the differential elliptic re- 
laxation equation for the velocity- pressure gradient contri- 
bution to fij. This can be obtained from the integral ex- 
pression (e.g. Manceau and Hanjalic, 2000) 


^j( x ) 

where 
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'M x ’ x ') = -«i( x )v 2 ^(x') ~ w j( x ) v 2 J^-( x ') 

+^S ij u k (x)V 7 -^-(x'), (12) 

and Gn(x,x') is approximated by the free-space Greens 
function = (47rr) -1 with r — ||x' — x||. The 

tensor function (s) (and <j){j ) can be represented by a 
polynomial expansion of basis tensors just as the associated 
Reynolds stress anisotropy tensor 6,j. For such a basis given 
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by T/ m - (x) (m = 1, . . . , n), the representation 


is assumed. A tensor scalar product (denoted by [ : ]) be- 
tween each basis tensor and the representation given in Eq. 
(13) can be formed, and this leads (using matrix notation 
for convenience) to 

[<A*(x):T (m) (x)] = j d 3 x' [**(x,x'):T (m >(x)] G n (x,x') 

= f d 3 x , 7n(x,x')7' nm Gii(x 1 x') (14) 

J n 


where 


T nm = [T(">(x):T( m >(x) 


(15) 


Since the functional dependency of [Tb l )(x):T^ m ^(x)] is 
only on x, Eq. (14) can be rewritten as 


[<£*( x ):T (m) (x) Tm = d 3 x'~/n (x, x')Gn( x , X T (16) 


where the scalar function 7 m (x,x') is modeled as 


7m ( x , x') = 7 m (x',x')exp 

If a tensor representation of the form 

N 

fw = E^ TW 

H = 1 



(17) 


(18) 


is assumed for </>*(x), then Eq. (16) can be rewritten as 

x f j3 , , , „exp(-r/L m ) 

Pm (x) — j a x 7 m (x , x ) — (19) 

This leads directly to the differential formulation 

(1 - L ly 2 ) /?m( x ) = (x, x) = 0 h m (x) (20) 

where (3^ (x) are the expansion coefficients from the ten- 
sor representation of the quasi-homogeneous form of the 
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pressure-strain rate correlation 0* 
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( 21 ) 


In this study, the pressure-strain rate model used in the el- 
liptic relaxation of the expansion coefficients differs slightly 
from the form used in the full differential form and given in 

Eq. (10) 


H h = - (c? + C\ j)e c b + KC 2 S 

+KC 3 (bS + Sb- |{bS}l) - KC 4 (bW - Wb) . 

( 22 ) 


A comparison of Eqs. (10) and (22) shows that the (slow) 
term proportional to b has been modified so that e c now 
affects the entire term. 

One of the improvements in the current elliptic re- 
laxation formulation is that the scaled velocity pressure- 
gradient correlation is 0(1) in the log-layer region. This 
rescaling negated the adverse influence of the elliptic oper- 
ator in the log-layer in the Durbin formulation. In order 
to insure the same benign effect of the elliptic operator on 
the expansion coefficients 0 rri} it is necessary to introduce a 
normalized set of basis tensors. 

In previous representations of the Reynolds stress an- 
isotropy tensor, the basis tensors used in a three-term ex- 
pansion were S, SW — WS, and S 2 — { S 2 } 1/3. In the 
log-layer, where the velocity gradient has a y~ l behavior, 
this then requires that the corresponding expansion coef- 
ficients 0 \ , and fis to have a y, y 2 , and y 2 behavior, 
respectively, in order to insure that </>* behaved as 0(1). 
Unfortunately, such a behavior would once again introduce 
an amplification effect by the elliptic operator on the ex- 
pansion coefficients. This problem can be easily remedied 
by introducing a scaled basis set of the form 

T (1) = S*, T (2) = S*W* - W*S*, 

T (3) = S* 2 - 1{S* 2 }I, (23) 

o 


where S* = S/jS 2 } 1 / 2 and W* = W/jS 2 } 1 / 2 . This scaled 
now' makes the behavior of both the expansion coefficients 
and basis tensors of 0(1) in the log-layer, wdiich then pre- 
cludes any adverse effect of the elliptic operator in the re- 
laxation equation given in Eq. (20). 

These results can be extended to the elliptic relaxation 
function f XJ . In the full differential form, the relaxation 
function fij is composed of both the velocity-pressure gra- 
dient component as well as a contribution from the Reynolds 


stress and dissipation rate anisotropies (see Eq. (2)). The 
corresponding elliptic relaxation equation for the expansion 
coefficients is then given by 


(1 - Lly 2 )!3 m = T ([n A :T< m) ] +2e c [b:T (m) ]) 
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where, for the most general formulation, each coefficient can 
have an associated length scale 
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The 0 m boundary conditions are listed in Table 2 as func- 
tions of fa. 

Table 2. 0 m BOUNDARY CONDITIONS. 
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DISCUSSION 


Solution & Calibration 

For the channel flow' used in this study the equations 
w’ere non-dimensionalized by the bulk viscosity and friction 
velocity at the wall (i.e. wall-units). The physical aspect of 
the channel flow was defined by the Reynolds number based 
on channel half-height, Re T . A one-dimensional finite- 
difference code was used for all computations. The differ- 
encing template was node-centered with clustering close to 
the wall using an exponential stretching function. The 500 
node grid had the first point at a height of 0.1 wall units. 
The channel Reynolds number determined the grid height. 
The K and e equations were implicitly coupled as were the 
T ij — fij equations. The variables U, K,e and the r,j were 
solved in a time-dependent, mode and were updated by the 
time step At. Time stepping was not used in the solution 
of the fij equations. 
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Comparison of the Models 

Figures 1 to 6 show the predictive performance of both 
the Tij - fij and t^j — f3 n formulations. As the figures show, 
the two methods are coincident with the only difference be- 
ing in the 22-component of the pressure strain rate correla- 
tion. (In this figure the DNS data is the unaltered ( r c — r) 
high Reynolds number results so that in the vicinity of the 
wall the data differs significantly from the computed re- 
sults.) This difference is due to the slight change in the 
defining relation for H h (cf. Eqs. (10) and (22)) used in 
the Tij — (3 n formulation. Even with this slight change, it 
is interesting to note that the placement of the composite 
dissipation rate only has a localized effect on the pressure- 
strain rate component. There were no noticeable differences 
in the other computed quantities such as the mean flow, 
turbulent kinetic energy or Reynolds stresses due to the po- 
sition of this term. 



y 

Figure 1. MEAN VELOCITY DISTRIBUTIONS. 



y 

Figure 2. SHEAR STRESS DISTRIBUTIONS. 
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Figure 3. TURBULENT KINETIC ENERGY DISTRIBUTIONS. 



y 


Figure 4. m NORMAL STRESS DISTRIBUTIONS. 


The consistency between the two formulations across 
the channel can be shown analytically. Since, for two- 
dimensional flows, a three-term basis set is sufficient to rep- 
resent a tensor function that only depends on the tensors S 
and W. This implies that the effect of the relaxation term 
sK f should be reproduced by a properly formulated relax- 
ation procedure for the expansion coefficients. Consider, 
for example, the right hand side of the f3\ equation (24) 
expanded and expressed in terms of the invariants formed 
from b, S and W 

([n' , :T (1) ] + 2 e c [brT^]) = - (c°+Clj^j £ c {bS'} 

+/vC 2 {SS*} + 2A'C7 3 {bSS*} 

— 2A'C 4 {bWS*} + 2e c {bS*}. (26) 
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Figure 5. r 22 NORMAL STRESS DISTRIBUTIONS. 



Figure 6. /7 22 COMPONENT OF THE PRESSURE-STRAIN RATE COR- 

RELATION. 

The only non-zero components of S and W yield 5* 2 = 
l/\/2 = IV* 2 . If these values are used in Eq. (26), the re- 
sulting expression is proportional to the corresponding be- 
havior of the right hand side of Eq. (6) for the f \2 compo- 
nent (/i 2 = /?iT (1) ), 

([n k :T (1 >] + 2e c [biT' 1 *]) = V2 (77f 2 + 2e c 6 12 ) (27) 

Since /i 2 = fii/\/2, Eq. (27) shows that the results from the 
fij formulation and the expansion coefficient formulation 
are equivalent. Similar results would be obtained for the 
other fin~ fij sets of equations. 

SUMMARY 

A methodology has been developed that introduces a 
polynomial representation of the tensor redistribution func- 


tion fij. An elliptic relaxation equation corresponding th) 
fij relaxation equation is formulated that allows ter'me 
relaxation of the polynomial expansion coefficients. The 
method is demonstrated fiere on a fully developed channel 
flow problem and the ^w^ nethods are found to give similar 
results. This consistency is further confirmed through an 
analysis of the right hand sides of the corresponding elliptic 
relaxation equations. Ths formulation is now being used in 
the development of a fully explicit algebraic Reynolds stress 
model that incorporates an elliptic relaxation equation. 
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